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Abstract. We discuss a novel approach to the mathematical analysis of equations with 
memory, based on the notion of a state. This is the initial configuration of the system at 
time t = which can be unambiguously determined by the knowledge of the dynamics 
for positive times. As a model, for a nonincreasing convex function G : IR + — > IR + such 
that 



G(0) 



lim G(s) > lim G(s) > 

s — >Q s— >oo 



we consider an abstract version of the evolution equation 

POO 

d tt u(x, t) - A G(0)u(x, t) + G'(s)u(x, t - s)ds 
L Jo 

arising from linear viscoclasticity. 
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1. Preamble 

1.1. A general introduction to equations with memory. Many interesting physical 
phenomena (such as viscoelasticity, population dynamics or heat flow in real conductors, 
to name some) are modelled by differential equations which are influenced by the past 
values of one or more variables in play: the so-called equations with memory. The main 
problem in the analysis of equations of this kind lies in their nonlocal character, due to the 
presence of the memory term (in general, the time convolution of the unknown function 
against a suitable memory kernel). Loosely speaking, an evolution equation with memory 
has the following formal structure: 

(1.1) d t w(t) = F(w(t),w t (-)), t>0, 

where 

w t (s) = w(t - s), s > 0, 

and T is some operator acting on w(t), as well as on the past values of w up to the actual 
time t. The function w is supposed to be known for all t < 0, where it need not solve the 
differential equation. Accordingly, the initial datum has the form 

w(t)=w (t), t<0, 

where w is a given function defined on (— oo, 0]. 

A way to circumvent the intrinsic difficulties posed by the problem is to (try to) rephrase 
(1.1) as an ordinary differential equation in some abstract space, by introducing an aux- 
iliary variable accounting for the past history of w, in order to be in a position to exploit 
the powerful machinery of the theory of dynamical systems. This strategy was devised 
by CM. Dafermos [9], who, in the context of linear viscoelasticity, proposed to view w l 
as an additional variable ruled by its own differential equation, so translating (1.1) into a 
differential system acting on an extended space accounting for the memory component. 

However, when dealing with (1.1), what one can actually "measure" is the function 
w(t) for t > 0. The practical consequences are of some relevance, since for a concrete 
realization of (1.1) arising from a specific physical model, the problem of assigning the 
initial conditions is not only of theoretical nature. In particular, it might happen that 
two different initial past histories wq lead to the same wit) for t > 0. From the viewpoint 
of the dynamics, such two different initial past histories are in a fact indistinguishable. 
This observation suggests that, rather than the past history w l , one should employ an 
alternative variable to describe the initial state of the system, satisfying the following 
natural minimality property: 

two different initial states produce different evolutions w(t) for t > 0. 

From the philosophical side, this means that the knowledge of w(t) for all t > determines 
in a unique way the initial state of the problem, the only object that really influences the 
future dynamics. 

Of course, the main task is then to determine, if possible, what is a minimal state 
associated to (1.1). Unfortunately, a universal strategy is out of reach, and the correct 
choice depends on the particular concrete realization of (1.1). Nonetheless, for a large 
class of equations with memory, where the memory contribution enters in the form of 
a convolution integral with a nonincreasing positive kernel, a general scheme seems to 
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be applicable. In this paper, we discuss an abstract evolution equation with memory 
arising from linear viscoelasticity, presenting an approach which can be easily extended 
and adapted to many other differential models containing memory terms. 

1.2. Plan of the paper. The goal of the next Section 2, of more physical flavor, is 
twofold. First, we present an overview on materials with hereditary memory, dwelling 
on the attempts made through the years to construct mathematical models accounting 
for memory effects. Next, we illustrate the physical motivations leading to the concept 
of minimal state representation. In Section 3, we introduce an abstract linear evolu- 
tion equation with memory in convolution form, which, besides its remarkable intrinsic 
interest, will serve as a prototype to develop the new approach highlighted in the pre- 
vious sections. After some notation and preliminary assumptions (Section 4), we recall 
the history approach devised by Dafermos (Section 5), whereas, in Section 6, we make 
a heuristic derivation of the state framework, which will be given a suitable functional 
formulation in the subsequent Section 7 and Section 8. There, we prove the existence of 
a contraction semigroup, whose exponential stability is established in Section 9, within 
standard assumptions on the memory kernel. In Section 10 and Section 11, we discuss 
the link between the original equation and its translated version in the state framework. 
Finally, in Section 12, we compare the history and the state formulations, showing the 
advantages of the latter. 

2. A Physical Introduction 

2.1. The legacy of Boltzmann and Volterra. The problem of the correct modelling 
of materials with memory has always represented a major challenge to mathematicians. 
The origins of modern viscoelasticity and, more generally, of the so-called hereditary 
systems traditionally trace back to the works of Ludwig Boltzmann and Vito Volterra 
[1, 2, 48, 49], who first introduced the notion of memory in connection with the analysis 
of elastic materials. The key assumption in the hereditary theory of elasticity can be 
stated in the following way. 

For an elastic body occupying a certain region B C M. N at rest, the deformation of the 
mechanical system at any point x G B is a function both of the instantaneous stress and 
of all the past stresses at x. 

In other words, calling u = u(x, t) the displacement vector at the point x E B at time 
t > 0, the infinitesimal strain tensor 



1 

£ — - 

2 



Vm + Vu 



T 



obeys a constitutive relation of the form 

e(x,t) = i(cr(x,t),a t (x,-)), 

where 

<t(x, t) and <t*(cc, s) = cr(x, t — s), 

with s > 0, are the stress tensor and its past history at (x,t), respectively. In the same 
fashion, the inverse relation can be considered; namely, 

cr(x, t) = or(e(aj, t), e t (x, •)). 
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The above representations allow the appearance of discontinuities at time t; for instance, 
e(x,t) may differ from \im s _, e t (x, s). In presence of an external force / = f(x,t), the 
related motion equation is given by 

d tt u(x, t) = V • cr(x, t) + f(x, t). 

The concept of heredity was proposed by Boltzmann, essentially in the same form later 
developed by Volterra within a rigorous functional setting. However, a more careful anal- 
ysis tells some differences between the two approaches. Quoting [39], "...when speaking of 
Boltzmann and Volterra, we are facing two different scientific conceptions springing from 
two different traditions of classical mathematical physics" . 

Boltzmann's formulation is focused on hereditary elasticity, requiring a fading initial 
strain history £°(x, •) = £ t (x, -)\t=o for every x G £>, i.e., 

(2.1) lim £°(x, s) = lim e(x, -s) = 0, 
so that, for every fixed (x,t), 

e(x,t) = I d£(x,y), 

J — oo 

and assuming the linear stress-strain constitutive relation at (x, t) in the Riemann-Stieltjes 
integral form 

(2.2) <r(x,t)= [ G(x,t-y)d£(x,y), 

J —CO 

where 

G = G(x,s), s>0, 

is a fourth order symmetric tensor (for every fixed s), nowadays called Boltzmann function. 
In particular, Boltzmann emphasized a peculiar behavior of viscoelastic solid materials, 
named relaxation property: if the solid is held at a constant strain [stress] starting from 
a given time to > 0, the stress [strain] tends (as t — > oo) to a constant value which is 
"proportional" to the applied constant strain [stress]. Indeed, if 

(2.3) e(x,t) = e(x,t ) = e (x), Vt > t , 
it follows that 

/to 
G(x, t — y)de(x, y) = G oo (x)£ (x), 
-oo 

where the relaxation modulus 

G cx> (x) = lim G(x,s) 

s— >oo 

is assumed to be positive definite. 

Conversely, the general theory devised by Volterra to describe the constitutive stress- 
strain relation is based on the Lebesgue representation of linear functionals in the history 
space. In this framework, he stated the fundamental postulates of the elastic hereditary 
action: 

• the principle of invariability of the heredity, 

• the principle of the closed cycle. 
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In its simpler linear version, the Volterra stress-strain constitutive relation reads 

(2.5) <r(x,t) =G (x)e(x,t)+ I G'(x, t - y)e(x, y)dy, 

where 

G {x) = limG(x,s) 

s— >0 

and the relaxation function G'(x, s) is the derivative with respect to s of the Boltzmann 
function G(x, s). It is apparent that (2.5) can be formally obtained from (2.2) by means 
of an integration by parts, provided that (2.1) holds true. In which case, the Boltzmann 
and the Volterra constitutive relations are equivalent. It is also worth noting that if (2.3) 
is satisfied for some t > and 

/to 
G'(x,t-y)e(x,y)dy = 0, 
-oo 

then, in light of (2.5), we recover the stress relaxation property (2.4), as 

lim cr(x, t) = G (x)e (x) + lim / G'(x, t — y)e (x)dy = G oo (x)£ (x). 

The longterm memory appearing in (2.2) and (2.5) raised some criticism in the scientific 
community from the very beginning, due to the conceptual difficulty to accept the idea 
of a past history defined on an infinite time interval (when even the age of the universe 
is finite!). Aiming to overcome such a philosophical objection, Volterra circumvented the 
problem in a simple and direct way, assuming that the past history vanishes before some 
time t c < (say, the creation time). Hence, (2.5) is replaced by 

(r(x,t) = G (x)e(x,t) + / G'(x, t - y)e(x, y)dy, 

Jt c 

and the motion equation becomes the well-known Volterra integro- differential equation 
d tt u(x,t) = V ■G (x)Wu(x,t) + W ■ [ G'(x,t-y)Vu{x,y)dy + f{x,t), 

Jt c 

which (besides appropriate boundary conditions) requires only the knowledge of the initial 
data u(x,t c ) and dtu(x,t c ). 

Remark 2.1. Here, we exploited the equality 

FS = ¥S T , 

which holds for any fourth order symmetric tensor F and any second order tensor S. 

2.2. Further developments: the fading memory principle. In the thirties, Graffi 
[30, 31] applied Volterra's theory to electromagnetic materials with memory, successfully 
explaining certain nonlinear wave propagation phenomena occurring in the ionosphere. 
Nonetheless, the modern theory of materials with memory was developed after World 
War II, when the discovery of new materials (e.g., viscoelastic polymers) gave a boost to 
experimental and theoretical researches. In the sixties, a lot of seminal papers appeared 
in the literature, dealing with both linear and nonlinear viscoelasticity [3, 8, 10, 36, 
37, 40, 41]. In particular, the thermodynamics of materials with memory provided an 
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interesting new field of investigations, mainly because some thermodynamic potentials, 
such as entropy and free energy, are not unique, even up to an additive constant, and 
their definition heavily depends on the choice of the history space (see, for instance, [5]). 

Along the same years, Coleman and Mizel [6, 7] introduced a main novelty: the notion of 
fading memory. Precisely, they considered the Volterra constitutive stress-strain relation 

(2.5) , with the further assumption that the values of the deformation history in the far past 
produce negligible effects on the value of the present stress. In other words, the memory 
of the material is fading in time. Incidentally, this also gave the ultimate answer to the 
philosophical question of a memory of infinite duration. The fading memory principle is 
mathematically stated by endowing the space (8 of initial strain histories e°(x,-) with a 
weighted L 2 -norm 

h{s)\e°{x,s)\ 2 dsdx, 

where the influence function h is positive, monotone decreasing, and controls the relax- 
ation function G' in the following sense: 

/i" 1 (s)|G / (a;,s)| 2 dsda; < oo. 

The theory of Coleman and Mizel encouraged many other relevant contributions in the 
field, and was the starting point of several improvements in viscoelasticity (see [11, 24, 29, 
33, 47] and references therein). On the other hand, as pointed out in [25, 26], the fading 
memory principle turns out to be unable to ensure the well-posedness of the full motion 
equation of linear viscoelasticity 

POD 

(2.6) d tt u(x,t) = V • G (x)Vu(x,t) + V • / G'(x, s)Vu(x,t - s)ds + f(x,t), 

Jo 

with known initial data u(x,0), dtu(x,0) and u(x,—s), for s > 0. Moreover, the arbi- 
trariness of the influence function h (for a given G') reflects into the non uniqueness of 
the history space norm topology. 

In order to bypass these difficulties, Fabrizio and coauthors [12, 16, 23, 24] moved in two 
directions. Firstly, they looked for more natural conditions on the relaxation function G', 
focusing on the restriction imposed by the second law of thermodynamics. Secondly, they 
tried to construct an intrinsically defined normed history space. To this end, in the spirit 
of Graffi's work [32], they suggested that any free energy functional endows the history 
space with a natural norm [17, 18]. In this direction, starting from [3, 10], many other 
papers proposed new analytic expressions of the maximum and minimum free energies 
[12, 15, 19, 20, 27]. The first and well-known expression of the Helmholtz potential in 
linear viscoelasticity is the so called Graffi- Volterra free energy density 

1 1 f°° 

* G (x,t) = -G^aOkOM)] 2 -- J (G\x,s)[e(x,t) - e\x,s)le(x,t) - e\x,s))ds, 

where G' is negative definite with s-derivative G" positive semidefinite. With this choice 
of the energy, the asymptotic (exponential) stability of the dynamical problem (2.6) with 
/ = has been proved, under the further assumption that, for some 5 > 0, the fourth 
order symmetric tensor 

G"(x,s)+SG'(x,s) 
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is positive semidefinite for (almost) every s > (see [9, 21, 28, 42, 43, 45]). Besides, the 
form of ^ G suggested the introduction of the new displacement history variable (see [9]) 

rf(x, s) = u(x, t) — u(x, t — s), s > 0, 

so that, with reference to (2.5), 

POO 

cr(x,t) = G OQ (x)e(x,t) - / G'(x, s)Vt7*(:e, s)ds. 

Jo 

Accordingly, the dynamical problem (2.6) translates into the system 

POO 

d tt u(x, t) = V • Goo(a;)Vw(a3, t) - V • / G'(x, s)Vrf{x, s)ds + f(x, t), 

d t rf(x, s) = d t u(x,t) - d s rf(x, s), 

which requires the knowledge of the initial data u(x,0), d t u(x,0) and rj°(x,s), where 
the initial past history rj°(x, s) is taken in the space f), dictated by ^/ G , of all functions 
77 = ij(x, s) such that 

n poo 

\\'n\\% = — I / (G\x, s)Vrj(x, s),Vr](x, s))dsdx < oo. 
Jb Jo 

2.3. The concept of state. Unfortunately, a new difficulty arises in connection with the 
above energetic approach. Indeed, depending on the form of G', different (with respect to 
the almost everywhere equivalence relation) initial past histories t]i,r]2 £ $) could produce 
the same solution to the motion problem (2.6) (clearly, with the same initial data u(x, 0) 
and d t u(x,0)). This is the case when 

n poo 

(2.7) / / G'(a;, s + t)V [rj^x, s) - r] 2 (x,s)]dsdx = 0, Vr > 0. 

Jb Jo 

Remark 2.2. As a consequence, there is no way to reconstruct the initial past history 
r)°(x, s) of a given material considered at the initial time t — 0, neither from the knowledge 
of the actual state of the system, nor assuming to know in advance the future dynamics. 

Noll [44] tried to solve the problem by collecting all equivalent histories, in the sense of 
(2.7), into the same equivalence class, named state of the material with memory. Never- 
theless, any two different histories in the same equivalence class satisfy the relation 

which implies that the history space 9) is not a state space (unless each equivalence class 
is a singleton) and is not a state function, as required by thermodynamics. Further 
efforts have been made to endow the state space of materials with memory with a suitable 
"quotient" topology, which is typically generated by an uncountable family of seminorms 
[13, 34]. In this direction, however, there is no hope to recover a natural norm on the 
state space. Indeed, the main obstacle consists in handling a space where each element is 
a set containing an infinite number of functions (histories). 

A different and more fruitful line of investigations was devised in [14] (see also [27]), 
through the introduction of the notion of a minimal state. Drawing the inspiration from 
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the equivalence relation (2.7), the authors called minimal state of the system at time t 
the function of the variable r > 

poo poo 

C\x, t) = - / G'(x, r + s) Vrf{x, s)ds = - G'(x, t + s) [e(x, t) - e(x, t - s)] ds. 
Jo Jo 

With this position, the stress-strain relation takes the compact form 

(2.8) <t(x, t) = G 00 (x)e(x, t) + C(x, 0). 

The difficulties of the previous approaches are circumvented: the minimal state space is 
a function space endowed with a natural weighted L 2 -norm arising from the free energy 
functional 

^{x,t) = -G^xMx,^ 2 - - J {{G')-\x,T)d T C\x,r),d T C\x,r))dT, 
which involves the minimal state representation. 

Remark 2.3. As a matter of fact, the history and the state frameworks are comparable, 
and the latter is more general (see Section 12 for details). In particular, as devised in [22] 
(see also Lemma 12.1), it can be shown that 

2.4. The problem of initial conditions. The classical approach to problems with mem- 
ory requires the knowledge of the past history of u at time t — 0, playing the role of an 
initial datum of the problem. This raises a strong theoretical objection: as mentioned in 
Remark 2.2, it is physically impossible to establish the past history of u up to time — oo 
from measurements of the material at the actual time, or even assuming the dynamics 
known for all t > 0. On the other hand, in the state formulation one needs to know the 
initial state function 

C (*,r) = C t (a!,r)| t= o, 

which, as we will see, is the same as knowing the answer of the stress subject to a constant 
process in the time interval (0, oo); namely, the answer in the future. At first glance, this 
appears even more conceptually ambiguous and technically difficult than recovering the 
past history of u. We will show that it is not so. To this end, let us rewrite equation (2.6) 
in the form 

(2.9) d tt u(x, t) = V- G (x)Vu(x,t) + V- [ G'(x,s)Vu(x,t-s)ds-F (x,t) + f(x,t), 

Jo 

having set 

POO 

F {x,t) = -V- G'(x,t + s)X7u(x,-s)ds 
Jo 

= V • [G(au, t)Vu(a;, 0) - G^x^uix, 0) - C°(au, *)] • 

Under, say, Dirichlet boundary conditions, and having a given assignment of the initial 
values u(x,0) and d t u(x,0), the problem is well-posed, whenever F is available (i.e., 
whenever £° is available). 
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Remark 2.4. The function F is not affected by the choice of the initial data, nor by 
the presence of the forcing term /. Moreover, if the initial past history of u is given, 
the above relation allows us to reconstruct F . In this respect, the picture is at least not 
worse than before. 

Remark 2.5. It is important to point out that, whereas in the previous approach the 
whole past history of u is required, here, in order to solve the equation up to any given 
time T > 0, the values of F (x,t) are needed only for t G [0,T]. 

Agreed that the hypothesis of an assigned initial past history of u is inconsistent, we 
describe an operative method to construct the function F by means of direct measure- 
ments, moving from the observations that materials with memory (such as polymers) are 
built by means of specific industrial procedures. Assume that a given material, after its 
artificial generation at time t — 0, undergoes a process in such a way that 



In which case, the equality 
holds. Indeed, 



u(x,t) = u(x,0), Vt > 0. 
C t (x,0) = C°(x,t) 



poo 

C f (x, 0) = - / G'(x, s) [e(x, 0) - e(x, t - s)] ds 
Jo 

G'(aj, s) [e(x, 0) - e(x, t - s)] ds 

pOO 

= - G'(x,t + s)[e(x,0) - e(x,-s)]ds = C°(x,t). 
Jo 

Thus, (2.8) entails the relation 

?(x,t) = (r(x,t)-G 0O (x)e(x,0), 

meaning that £°, and in turn F , can be obtained by measuring the stress <r(x,t), for all 
times t > 0, of a process frozen at the displacement field u(x, 0). 

Next, we consider the equation of motion (2.9), relative to a material generated by 
the same procedure, but delayed of a time > 0. For this equation, the corresponding 
function F (x,t) is now available. 

Remark 2.6. As a matter of fact, F (x,t) is not simultaneously available for all t > 0. 
However, since the first process (which constructs F ) keeps going, at any given time 
T > 0, referred to the initial time t = of the problem under consideration, we have the 
explicit expression of F (x,t) for all t G [0,£d + T], which is even more than needed to 
solve (2.9) on [0,T]. 

3. An Abstract Equation with Memory 

We now turn to the mathematical aspects of the problem, developing the state approach 
for an abstract model equation. 

Let H be a separable real Hilbert space, and let A be a selfadjoint strictly positive 
linear operator on H with compact inverse, defined on a dense domain T>(A) C H . For 
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t > 0, we consider the abstract homogeneous linear differential equation with memory of 
the second order in time 

rt 



(3.1) d tt u(t) + A au(t) - / fi(s)u(t - s)ds 







= 0. 



Here, a > 0, £ G (0, oo] and the memory kernel 

n : n = (o,e) -> (o,oo) 

is a (strictly positive) nonincreasing summable function of total mass 

ri 

[i(s)ds E (0, a), 



I 

Jo 



'0 

satisfying the condition (automatically fulfilled if £ = oo) 

lim /i(s) = 0. 

The dissipativity of the system is entirely contained in the convolution term, which ac- 
counts for the delay effects: precisely, finite delay if £ < oo, infinite delay if £ = oo. The 
equation is supplemented with the initial conditions given at the time t — 

(u(0) = m , 
d t u(0) = v , 
u{s)\sen = 0o(s), 
where uo, fo and the function <fi , defined on Q, are prescribed data. 

Remark 3.1. A concrete realization of the abstract equation (3.1) is obtained by setting 
= IR + , H = [L 2 (B)] N , where B C M. N is a bounded domain with sufficiently smooth 
boundary dB, and 

A = -A with V(A) = [H 2 (B)] N n [H^(B)] N . 

In that case, calling 

a — fx(a)da = G(s) 
Jo 

and 

u(t) = u(x, t), x E B, 

the equation reads 

d tt u(x, t) - A G(0)t*(au, t) + J G'(s)u(x, t - s)ds] = 0, 
u(x,t)\ xedB = 0, 

and rules the evolution of the relative displacement field u in a homogeneous isotropic 
linearly viscoelastic solid occupying a volume B at rest [24, 47]. 

Putting /i(s) = if s > £, and defining 

(3.3) F (t) = [ n{t + s)0 o (s)ds, 

Jo 
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equation (3.1) takes the form 

(3.4) d tt u{t) + A au(t) - [ n(s)u(t - s)ds - F (t) =0. 

L Jo - 1 

Introducing the Hilbert space 

V = V(A^ 2 ), 
with the standard inner product and norm 

(ui, u 2 ) v = (A l/2 Ul , A 1/2 u 2 ) H , \\u\\ v = \\A l/2 u\\ H , 
we stipulate the following definition of (weak) solution. 

Definition 3.2. Let u G V, v G H and (f) : ^ V be such that the corresponding 
function F given by (3.3) fulfills 

F eLU[0,oo);V). 

A function 

u G C([0, oo),V)n C^fO, oo), H) 
is said to be a solution to the Cauchy problem (3.1)-(3.2) if 

u(0) = u , d t u(0) = v , 

and the equality 

(d tt u(t),w) + a(u(t),w) v - / pL(s){u(t - s), w) v ds - {F (t), w) v = 

Jo 

holds for every w eV and almost every t > 0, where (•, •) denotes duality. 

4. Notation and Assumptions 

4.1. Notation. The symbols (•, -)x and || • ||x stand for the inner product and the norm 
on a generic Hilbert space X, respectively. In particular, for the spaces H and V, we have 
the well-known norm relations 

\\u\\ v = \\A 1/2 u\\ H > ^/\~ 1 \\u\\ H , \/u G V, 

where Ai > is the first eigenvalue of A. We denote by 

V* = V(A~ 1/2 ) 

the dual space of V, and by (■, •) the duality product between V* and V. We also recall 
the equality 

\\ w \\ v . = \\A-^ 2 w\\ H , Vw G V*. 

For a nonnegative (measurable) function u on fi = (0,^) and for p = 1,2, we define the 
weighted L p -space of X-valued functions 

ft 

Ll{n-X) = [ij:n^X : J u(s)U(s)\\ p x ds < oo}, 

normed by 

; ( n^) = (j^a;( a ) ll^l^da) " 
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If p = 2, this is a Hilbert space endowed with the inner product 

(V'1,'02)l2(Q;X) = / Uj{s){lp 1 (s),1p 2 (s))xds. 

Jo 

Finally, given a generic function ip : Q — > X, we denote by its distributional derivative. 

.A word of warning. In order to simplify the notation, if tp is any function on Q, we agree 
to interpret if)(s) = whenever s £ Q (in particular, if £ < oo, whenever s > £). 

4.2. Assumptions on the memory kernel. As anticipated above, /i : — > (0, oo) is 

nonincreasing and summable. Setting 

M(s) = J n(cr)dcr = J fj,(s + a)da, 
we require that M(0) < a. 

Remark 4.1. Note that M(s) > for every s £ [0,£), and lim s ^M(s) = 0. 

For simplicity, we will take 
(4.1) a-M(0) = l. 

In addition, we suppose that /x is absolutely continuous on every closed interval contained 
in Q. In particular, fi is differentiable almost everywhere in Q and // < 0. Finally, /i is 
assumed to be continuous at s — £, with /i(£) = 0, if £ < oo. Conversely, if £ = oo, as 
H is nonincreasing and summable, we automatically have that fj,(s) — > as s — > oo. In 
fact, we could consider more general kernels as well, allowing fi to have a finite or even 
a countable number of jumps (cf. [4, 45]). However, in this work, we will restrict to the 
continuous case, in order not to introduce further technical difficulties. 

5. The History Approach 

An alternative way to look at the equation is to work in the so-called history space 
framework, devised by Dafermos in his pioneering paper [9], by considering the history 
variable 

rf(s)=u(t)-u(t-s), t>0,sGfi, 
which, formally, fulfills the problem 

d t rf{s) = -d s rf(s) + d t u(t), 
t/(0) = 0, 
r)°(s) = u - (p (s). 

To set the idea in a precise context, let us introduce the history space 

M = Ll(Q;V), 

along with the strongly continuous semigroup R(t) of right translations on J\4, namely, 



(R(t)r])(s) 



< s < t, 

r)(s — t) s > t, 
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whose infinitesimal generator is the linear operator T defined as (cf. [35, 45]) 

T V = -Dr), V{T) = {rj G M : Dr] G M, r](0) = 0}, 

where 77(0) = lim s ^or)(s) in V. Then, recalling (4.1), equation (3.1) translates into the 
differential system in the two variables u = u(t) and r\ = r/(s) 



(5.1) 



fd tt u{t)+A u(t) + J fx(s)r ] t (s)ds 
[dtrf = Trf + d t u(t). 

Accordingly, the initial conditions (3.2) turn into 



0, 



{U(0) = M , 
d t u(0) = v , 
V° = Vo, 

where 

(5.3) r}o(s) = u - <j>o{s). 

Introducing the extended history space 

m = v x H x M, 

normed by 

IIK^)IIot= \W\\v+ \\ v \\h+ WvWm, 

problem (5.1)-(5.2) generates a contraction semigroup S(t) on 9JI (see [21, 35, 45]), such 
that, for every (u ,v ,r]Q) G 971, 

T,(t)(u 0: v 0: rj ) = {u{t),d t u{t),rf). 

Moreover, rf has the explicit representation 



(5.4) if{s) = 



u(t) -u(t-s) < s < t, 

Vo( s ~ + u (t) — u o s > t. 



Concerning the relation between (5.1)-(5.2) and the original problem (3.1)-(3.2), the fol- 
lowing result holds [35]. 

Proposition 5.1. Let (uo,Vo,i] ) G 971. Then, the first component u(t) of S(t)(-u , v , rj ) 
solves (3.1)-(3.2) with 



F (t)= [ fi(t + s){u -r] (s)}ds. 
Jo 

It is easy to see that i] G M implies that F G L°°(]R + ; V). 
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6. The State Approach 

An essential drawback of the history approach is that, for given initial data uq and t> , two 
different initial histories may lead to the same solution u(t), for t > 0. Somehow, this is 
not surprising, since what really enters in the definition of a solution to (3.1)-(3.2), rather 
than 0o (which, by (5.3), is related to the initial history 770), is the function F , defined 
in (3.3) and appearing in equation (3.4). Thus, from the dynamical viewpoint, two initial 
data (f) i and 002 should be considered by all means equivalent when the corresponding 
function Fqi and F02 coincide, due to the impossibility to distinguish their effects in the 
future. On this basis, it seems natural to devise a scheme where, rather than 0o, is the 
function F to appear as the actual initial datum accounting for the past history of u. 

In order to translate this insight into a consistent mathematical theory, it is quite 
helpful to see first what happens at a formal level. To this aim, for t > and r G (], we 
introduce the (minimal) state variable 

C'(r)= / v(t + s){u(t) - u(t - s)}ds, 
Jo 

which fulfills the problem 

'd t C t (r) = d T C(r) + M(r)d t u(t), 

CO = 0, 

.C°(r)=Co(r), 

having set 

Co(t)= / fi{T + s){u -Ms)}ds = M{T)u -F {r). 
Jo 

Accordingly, in light of (4.1), equation (3.1) takes the form 

duu^+Alu^ + CiO)} =0, 
where ^ 

C'(0) = limC*(r)= / »(s){u(t)-u(t-s)}ds. 
Jo 

Rather than £*, it seems more convenient to consider as a state the new variable 

e(r) = -d T C\r) = - f V(r + s){u(t) - u(t - s)}ds, 

Jo 

which, in turn, fulfills the problem 

'd t e(r)=d T e(r)+^r)d t u(t), 

where the initial datum £ reads 

6( T ) = - / /A r + s){uo - 0o(s)}ds = h{t)uq + / fx'(T + s)(j) (s)ds. 
Jo ' Jo 

If £ < 00, we have also the "boundary" condition 

= 0, 



A NEW APPROACH TO EQUATIONS WITH MEMORY 15 

which comes from the very definition of £*. Since (*(£) = 0, we find the relation 
(6.1) /V(r)dr = C'(t ), Vr G Q. 

" TO 

In particular, in the limit r — > 0, 

/V(r)dr = C*(0). 
Jo 

Therefore, (3.1)-(3.2) is (formally) translated into the system 



(6.2) 

with initial conditions 



{d tt u(t) + A[u(t)+ J e(r)dr 
l^(r) = W) + MrM), 



0, 



(«(0) = M , 
d t u(0) = v , 

Remark 6.1. Observe that the nonlocal character of (3.1) is not present in (6.2) any 
longer, since it is hidden in the new variable £*. 

At this point, to complete the project, two major issues need to be addressed: 

• Firstly, we have to write (6.2)-(6.3) as a differential equation in a suitable func- 
tional space, providing an existence and uniqueness result. 

• Secondly, we have to establish a correspondence (not only formal) between the 
solutions to (6.2)-(6.3) and the solutions to the original problem (3.1)-(3.2). 

7. The State Space 

The first step to set (6.2)-(6.3) in a proper functional framework is to interpret in a correct 
way the derivative d T appearing in the second equation of (6.2). We introduce the new 
memory kernel 

^( r ) = -L:Q^[0,oo), 
Mr) 

and we put 

z/(0) = lim v{t). 

In view of the assumptions on /i, the function v is continuous and nondecreasing on Q, 
with nonnegative derivative (defined a.e.) 

Ar)= 



WW 

Moreover, 

lim v(t) = oo. 

Introducing the state space 

V = Ll(n;V), 
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whose norm is related to the free energy functional \1/f of Section 2.3, we consider the 
strongly continuous semigroup L{t) of left translations on V, defined by 

(L(0O(t)=£(* + t). 

It is standard matter to verify that the infinitesimal generator of L{t) is the linear operator 
P on V with domain 

V{P) = {ieV :D(GV, m=0}, 
where £(£) = lim r _^£(r) in V, acting as 

Pt = Dt, ^eV(P). 

Note that, if £ G V(P), then \\^\\ v G C(Q). 

Remark 7.1. If £ = oo, the condition £(oo) = is automatically satisfied whenever 
£, D£ G V. Indeed, for every r > 1, the function £ is absolutely continuous on [l,r] with 
values in V, and 



£(r) = £(1) + y L>£(s)ds. 

On the other hand, 

f \\Dt{s)\\ v te = £ v 7 ^) \[W)\\Dt{s)\\v&s 

< (^ r M^)d S ) 1/2 (^ r K S )P£( S )||^d S ) 



1/2 



< ^/M(l)|p£|| v . 

lim / D£;(s)ds 



£(oo) = lim £(r). 



Thus, 

exists in V, and so does 
Since the function 

T » "{t)\\Z{t)\\$ 

is summable and u{t) — > oo as r — > oo, it must necessarily be £(oo) = 0. Arguing in a 
similar manner, we see that 

CZ^GV =► G L~(fi). 

Lemma 7.2. For every £ G T>(P), 

u\r)U(r)\\ 2 v dr<oc, 



and the limit 

exists finite (equal to zero if v(0) = 0). Moreover, 

(7-1) 2(P£,£) V = -i/(0)||e(0)||^ - /"VwiieWll^dr < 0. 

Jo 
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Proof. We begin to prove the existence of a sequence £ n j £ such that 

"(OIK(Ollv->o. 

If £ — oo, this is a direct consequence of the summability of ^||£||y. Conversely, if £ < oo, 
for every t < £ we have (recalling that = 0) 

"(r)||£(r)|&< ( jf ll^(*)llvd«) 2 < (€-r)||^||=. 

Let now e n < £ n be any sequence such that e n { 0. Then, 

2<P£,0v= lim r^(r)A|| e ( r) ||^dr 
n^oo / ar 

= hm { -u(e n ) W&nWv- t u'(r)U(r)\\ 2 v dr}- 

Since the limit exists finite, both summands have the same sign and the integral term is 
monotone, we conclude that the limit of the sum equals the sum of the limits, so yielding 
equality (7.1). We are left to demonstrate the implication 

1/(0) = =► u(0)\\m\\ 2 v = o. 
Indeed, choosing an arbitrary r G f2, for any r < r we have 

u(t) U(r) \\l < 2u(t) U(t ) \\ 2 y + 2i/(r) ( jH || v d S ) * 

< 2i/(r) ||e(ro) f y + 2 ( jf ^ \\ v ds) * 

<2u(r)\\ar )\\ 2 v + 2(r -r)\\Dal 

Therefore, 

]imsn V u(r)U(r)\\ 2 v <2r \\DC\\ 2 v , 
and letting r — > 0, the claim follows. □ 

The following simple lemma will be needed in the sequel. 
Lemma 7.3. Let (eV. Then, £ E L 1 ^; V) and 

' l \\i(r)\\ v dr<^MW)U\\v- 



As a byproduct, the function 

t ^ J £(r)dr 

belongs to C([0, oo), V) and vanishes at infinity. 
Proof. Using the Holder inequality, 

U(r)\\vdr= f e ^)V^)\\ar)\\vdr<VM(^U\\v, 



I 

Jo 



10 

as claimed. □ 
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8. The Semigroup in the Extended State Space 

We are now in a position to formulate (6.2)-(6.3) as an abstract evolution equation on a 
suitable Hilbert space. To this end, we introduce the extended state space 

21 = v x H x V, 

normed by 

IIK^OII^IHI 2 y + IHII + lffi, 

and the linear operator A on 21, with domain 

V(A) = {(u,v,Z)e%:veV,u + f^(r)dr G V(A), £ G P(P)}, 

acting as 

A(u, v, = (u, -A [u + J £ e(r)dr] , P£ + H • 
Introducing the 3-component vectors 

Z(*) = (u(i), u(f), £*) and z = (w , vq, f ) e 21, 
we view (6.2)-(6.3) as the Cauchy problem in 21 

(8.1) (^(«) = AZ( t ), 

The following result establishes the existence and uniqueness of a (mild) solution 

Z G C([0,oo),2l). 

Theorem 8.1. Problem (8.1) generates a contraction semigroup S(t) = e tA on 21 such 
that 

Z(t) = S{t)z, Vt > 0. 

Moreover, the energy equality 

(8-2) f t \\S(t)z\\l = -^(0)11^(0)11^ - jy'(r)\\e(r)\\ldr 

holds for every z G D(A). 

Proof. On account of the classical Lumer-Phillips theorem [46], we know that A is the 
infinitesimal generator of a contraction semigroup on 21 provided that 

(i) the inequality (Az, z)% < holds for every z G T>(A); and 

(ii) the map I — A : V(A) ^ 21 is onto. 
Concerning point (i), from (7.1) we see at once that 

(A*,z) a =(P£,Ov<0, 

for every z = (u,v,£) G P(A). 

In order to prove (ii), let z* = (u*,v+,£+) G 21 be given. We look for a solution 
z— (u,v,£) G V(A) to the equation 
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(8.3) 



u — V = u+, 



v + A 



u 



+ [ f (r)drj = v*, 
Jo J 



t(T)-Dt(T)-IM(T)v = MT)- 

Given a function g on VL (extended on the whole real line by setting g(s) = if s G" fl) 
and denoting 

= e s X(-oo,o](s), 

we consider the convolution product in K. of S and g at the point r G fl 



(S*g)(r)= S(r-s)g(s)ds= / e T -g(s)ds. 
Jr Jt 

It is well known (see, e.g., [38]) that 

geL 2 (fl) =>- £*geL 2 (fl) 

and 

\\£ * 9\\l*(si) < IMU 2 («)- 
For any fixed i> G V, we define the function 

(8.4) t(T)=v(E*IM)(T) + (e*S il )(T). 

We begin to show that £ G V. Indeed, 

r^)||e(r)|| 2 y dr<2|| W || 2 y f v{r)\{£^){r)\ 2 dr + 2 ( v{r) \\ (£ * &)(r) f v dr 
Jo Jo Jo 

' _ / \ 2 

e r s V z/ ( r ) A*(s)ds ) dr 

e T - s v^)||e,(s)||yd S ) 2 dr 
< 2||u||^||f * ||i 2(Q) + 2||£: * (v^ ||£*||v)||i2 (n) 



< 2|k||y 

+ 2 



<2M(0)||t;|| 2 / + 2||e,||^. 



Moreover, 

Since (cf. Lemma 7.3) 
we conclude that 



U(r)\\ v < I <?- s {»{s)\\v\\ v +\\Us)\\vMs. 

s^ii(s)\\v\\ v +\%{s)\\v eL\n), 
lim||e(r)|| v = 0. 



Taking the distributional derivative in both sides of (8.4), it is apparent that such a £ 
satisfies the third equation of (8.3). Moreover, by comparison, it is readily seen that 
-D£ G V. In summary, £ G V(P) (for any given v E V) and fulfills the third equation of 
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(8.3). At this point, we plug £ into the second equation of (8.3), reading u from the first 
one. Noting that 

7 = 1 + ^ (J e r ->(s)ds)dr = 1 + J fi(s){l - e' s }ds > 0, 

we obtain 

(8.5) v + jAv = v i ,-A(u i , + w), 

having set 

Jo 

The elliptic equation (8.5) admits a (unique) solution v G V, provided that its right-hand 
side belongs to V*, which immediately follows from w G V. Indeed, using the Holder 
inequality, 

\H\v< f \\(S*^)(r)\\vdr 
Jo 

< 



w 



[ (S * \\Uv)(r)dr 
Jo 

< [ ^){S*(V^\\U\v))(r)dr 
Jo 



< y/M(ti)M*\\v- 

Finally, by comparison, we learn that 



u+ f £(r)dr = - v) G V{A). 

Jo 



This completes the proof of point (ii). 

We now appeal to a general result of the theory of linear semigroups [46]. Namely, if 
z G £>(A), then 

S{t)zeV(A), Vt>0, 

and 

j t \\S(t)z\\l = 2(AS(t)z,S(t)z) sl . 

On the other hand, since 

(AS(t)z,S(t)z)a={PZ t ,Z t ) v , 
the energy equality (8.2) follows from Lemma 7.2. □ 

Corollary 8.2. The third component of the solution S(t)z (the state) has the explicit 
representation formula 

(8.6) £'(t) = fo(* + r ) + A*0>(*) - p(t + r)u + f /i'(r + s)u(t - s)ds, 

Jo 

which is valid for every z = (uq, vq, £o) £ 21- 
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Proof. Assume first that z lies in a more regular space, so that d t u e L\ oc ([Q,oo);V). 
Then, £' satisfies the nonhomogeneous Cauchy problem in V 

Applying the variation of constants to the semigroup L(t) = e tp (see [46]), we obtain 

£'(r) =&(t + r)+ f fi(t + t - s)d t u(s)ds. 
Jo 

The desired conclusion (8.6) is drawn integrating by parts. Using a standard approxima- 
tion argument, the representation formula holds for all z e 21. □ 

Remark 8.3. In the above corollary, the continuity properties of fx play a crucial role 
when integrations by parts occur. Nonetheless, if fx has jumps, it is still possible to find 
a representation formula, which contains extra terms accounting for the jumps of fx. 

Remark 8.4. The state variable £' is minimal in the following sense: if (u(t) , dtu(t) , £*) 
is a solution to (8.1) with -u(t) = for every t > 0, then is identically zero. Indeed, on 
account of (8.1) and (8.6), 

e(r)=Ut + r), Vt>0, 

and 

/■^ ft ft 

0= / ^(r)dr= / eo(t + r ) dr = / ^ ( r )dr, Vt > 0, 

which implies that £o = and, in turn, £* = 0. 

9. Exponential Stability 

9.1. Statement of the result. We prove the exponential stability of the semigroup S(t) 
on 21, assuming in addition that fx satisfies 

(9.1) fu'(s) + 5fi(s) < 0, 
for some 5 > and almost every sfO. 

Theorem 9.1. Let fx satisfy (9.1). Then, there exist K > 1 and u > snc/i £/ia£ 

(9.2) ||S(^|| a <X|H| a e-"*, 
/or every z G 21. 

Before proceeding to the proof, some comments are in order. Condition (9.1) is quite 
popular in the literature; indeed, it has been employed by several authors to prove the 
exponential decay of semigroups related to various equations with memory in the history 
space framework (e.g., in connection with the present equation, [21, 28, 42, 43]). On the 
other hand, the recent paper [45] shows that the exponential decay for such semigroups 
can be obtained under the weaker condition 

(9.3) fx(a + s) <Ce- 5rT fi(s), 
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for some C > 1, every a > and almost every sffl, assuming that the set where fi' = 
is not too large (in a suitable sense). It is apparent that (9.3) and (9.1) coincide if C = 1. 
However, if C > 1, then (9.3) is much more general. For instance, it is always satisfied 
when £ < oo (provided that ji fulfills the general assumptions of Section 4). On the 
contrary, (9.1) does not allow fi to have flat zones, or even horizontal inflection points. As 
shown in [4], condition (9.3) is actually necessary for the exponential decay in the history 
space framework. This is true also in the state framework. 

Proposition 9.2. Assume that the semigroup S(t) on 21 is exponentially stable. Then, fi 
fulfills (9.3). 

We omit the proof of the proposition, which can be obtained along the lines of [4], 
showing that the exponential stability of S(t) implies the exponential stability of the 
left-translation semigroup L(t) on V. 

We finally point out that, although we stated the theorem using (9.1), the result is still 
true under the more general hypotheses of [45] (but a much more complicated proof is 
needed) . 

9.2. Proof of Theorem 9.1. Appealing to the continuity of S(t), it is enough to prove 
inequality (9.2) for all z G V(A). Fix then 

z = (w ,^o,£o) e T>(A), 



For an arbitrary (3 G fi, we define the (absolutely continuous) function p : Q — > [0, 1] 



and denote 



S(t)z 



(u(t),v(t),e)eV(A). 



Introducing the energy 



E{t)= l -\\S(t)z 





and we consider the further functionals 




* 2 (t) 



{v(t),u(t)) H . 



Recalling Lemma 7.3, 



(9.5) 
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Thus, from the continuous embedding V C H, 
(9.6) \$ t (t)\<CoE(t), ^ = l,2, 

for some Cq > independent of the choice of z G X>(A). 
Lemma 9.3. There is c\ > independent of z such that 

(9-7) ^(t) < M(f3){±\\u(t)\\l - ±\\v(t)\\H + ciMYv}. 

Proof. We have 



d^ 1 



f p(r)(d t v,ar)) H dr- [ p(r)(v, d t ^r)) H dr. 
o Jo 



We now estimate the two terms of the right-hand side, exploiting the equations of (8.1) 
and the integral control (9.5). For the first one, 

ft pi ft ft 

' p(r){d t v,i(r)) H dr = J p(r)(u, £(r)) y dr + J p(r) ( J (£(/), £(r)> y dr')dT 

^ll^r)||ydr+(^||e(r)||ydr) 2 



< u 



< v/M(o)ii«i| v |ieiiv+M(o)iieii^ 

<^M(P)\\u\\l + M(0)[l + J ^)u\\ 2 v . 

Concerning the second term, we preliminarily observe that, since £ G T>(P) for all times, 
we have (cf. Remark 7.1) 

sup||£(r)|| v <oo and ||^)||y = 0. 
Thus, an integration by parts gives 

- J p(r) (v, P^r)) H dr = - J p(r) A (v, £(r)) H dr = i jf <„, £(r)>„dr. 



Hence, 



■y p(t)(»,^(t)Mt = -(| p(r)Mr)dr) 11* + - jf (^(r)) fl dr 

< -M^Ht;!^ + -A- Nk / lieWUvdr 



,2 , V^W 



<-M(/?)|M|^ + ^^|MM|£||y 



Collecting the above inequalities, the conclusion follows. □ 
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Lemma 9.4. The functional $2^) fulfills the differential inequality 

(9-8) ±* 2 (t) < -?||u(f)|& + \\v(t)\\ 2 H + M(0)\\e\\l 

Proof. By virtue of (8.1) and (9.5), 

d f e 

= ~\\u\\l + \\v\\% - (u^(r)) v dr 



<-\H\l + \M% + y/M^)\\u\\vU\\v 

<-lh\\ 2 v + \M 2 H + M(o)U\\ 2 v, 

as claimed. □ 
At this point, we define the functional 

•0 = ^.(0 + *.«). 

which, due (9.7) and (9.8), satisfies the differential inequality 

(9.9) + E(t) < ciUTv, 
for some c 2 > independent of z. Besides, in light of (9.6), 

(9.10) mt)\<c 3 E(t), 
with c 3 = c (3/M(/3) + 1). Finally, we fix 



. r S li 

£ = mm < - — , - — > 
I 2c, 2c, J 



and we set 

= £(t) + e$(t). 

Note that, by (9.10), 

and in turn, by (9.4) and (9.9), 

— tf(i) + 2wtf(t) < 0, 

with a; = e/3. Therefore, the standard Gronwall lemma yields 

\\S(t)z\\l = 2E(t) < 4tf(*) < 4^(0)e~ 2 ^ < 6£(0)e~ 2 ^ = 3 1 1 ^ 1 1 ^ e~ 2 ^* . 
The proof of Theorem 9.1 is completed. 

Remark 9.5. Observe that the proof of Theorem 9.1 is carried out employing only energy 
functionals, and it makes no use of linear semigroup techniques. Thus, the same energy 
functionals can be exploited to analyze semilinear versions of the problem (for instance, 
to prove the existence of absorbing sets and global attractors). 
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10. The Original Equation Revisited 

Somehow, this novel state approach urges us to consider the original problem under a 
different perspective. Indeed, as we saw in Section 6, the solutions to (3.1)-(3.2) are 
determined, besides by Uq and Vq, by the knowledge of the function F , and not by the 
particular form of the initial past history O - Therefore, with reference to Definition 3.2, 
we introduce the class of admissible past history functions 

A = {0 : n - V : t h- J\{t + S )0( S )ds G LU[0, 00); V)}, 
and we define the linear map 

A:A^LU[0,oo);V) 

as 

<f> h+ A<j>(t) = / n(t + s)(f)(s)ds. 
Jo 

Note that A0(t) = if t > t. Accordingly, we define the class of state functions 

S = AA. 

Clearly (and this is really the point), the map A may not be injective, meaning that 
different <fi G A may lead to the same element of S. 

Coming back to Definition 3.2, the assumption on F can now be rephrased as 

F = A0 O with 0o e A, 

and we can reformulate the definition of solution to (3.1) in the following more convenient 
(and certainly more physical) way. 

Definition 10.1. Let the triplet 

(^0,^0,^0) eV x H x S 

be given. A function 

u G C([0, 00), V) n C^QO, 00), H) 
is said to be a solution to equation (3.1) with initial state (uo,v ,F ) if 

u(0) = u , d t u(0) = v , 

and the equality 

(d tt u(t), w) + a(u(t), w) v - / fi(s)(u(t - s), w) v ds - (F (t), w) v = 

^0 

holds for every w G V and almost every t > 0. 

In this definition, the initial datum O has completely disappeared, since the state func- 
tion Fq contains all the necessary information on the past history of the variable u needed 
to capture the future dynamics of the equation. Hence, we removed the (unphysical) am- 
biguity caused by two different initial histories leading to the same state function, which, 
as we saw, is what really enters in the definition of a solution. 
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Remark 10.2. We point out that the function F (t) is not influenced by the dynamics 
for t > 0, nor by the presence of a possible external force. As a matter of fact, if the initial 
past history <J)q is known, then Fq is uniquely determined by (3.3). On the other hand, 
even if the particular O leading to F Q is unknown, in principle, F can still be determined 
(cf. Section 2.4). 

The remaining of the section is devoted to investigate the properties of the space S. 
We begin with a lemma, which provides a precise formulation of the formal equality (6.1), 
devised in Section 6. 

Lemma 10.3. Whenever (p G A, the map 

t» [ //(r + s)||0(s)|| y ds 

Jo 

belongs to L 1 (t, oo) for every t > 0, and the equality 

ft ft 

(10.1) k<j>(t) = -J (y /i'(r + S )0( S )ds)dr 

holds for every t > 0. Moreover, if <fi G L^Q; V), then <fi G A and (10.1) holds for every 
t > 0. 

Proof. Let G ^4 be given. For every fixed t > 0, 

A0(t o ) G V, for some t <t. 

Since /x is a nonincreasing function and A0(t o ) is a Bochner integral, this is the same as 
saying that 

/ fi{t + s)\\<f){s)\\ v ds < I fj,{t + s)\\<j)(s)\\ v ds < oo. 
Jo Jo 
Exploiting the equality 

fj,(t + s) — — J (j,'(t + s)dr, 

and exchanging the order of integration, we conclude that 
ft ft ft 

J ^{t + s)U(s)\\ v ds = - y (jf /i'(r + S )||0( S )|| y d S )dr<oo. 

Hence, 



/ /i'(r + s)||0(s)||yds G oo), 

and (10.1) follows from the Fubini theorem. Concerning the last assertion, just note that 
G LliVL; V) if and only if A0(O) G V. □ 

Remark 10.4. As a straightforward consequence of the lemma, 

S C C ([t,oo),V), Vt>0, 

where Co denotes the space of continuous functions vanishing at infinity. 

Given F G <S, it is then interesting to see what happens to F(t) in the limit t — > 0. 
Three mutually disjoint situations may occur: 
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(i) \im t ^ F(t) exists in V; 

(ii) F G L QO (R + ; V) but lim^o does not exist in V; 

(iii) ||F(t)||y is unbounded in a neighborhood of t — 0. 

As we will see, (i) is the most interesting case in view of our scopes. For this reason, we 
introduce the further space 

S = {F G S : 3 limF(t) in V\. 

1 t^o ' 

In light of Remark 10.4, it is apparent that 

<S cCo([0,oo),n 

We preliminary observe that if F = A0 with G V), then Lemma 10.3 yields at 

once 

limF(t) = A0(O) in V, 

so that F G «So- However, the picture can be more complicated. Indeed, it may happen 
that F G So but A0(O) is not defined for any G A _1 F, as the following example shows. 

Example 10.5. Consider the kernel 

H(s) = l-s, n=(o, 1). 
Given any nonzero vector u G V, set 



F(t) 



t) sin 1 



sin — dx 

x 



X[o,i](t) u, 



which clearly satisfies 
Then, F = A0 with 



limF(t) = 



= - 



sin 1 
1 



/ sin —dx 
Jo x ■ 



u. 



1 



LI 



cos 



1 - si 



u. 



but A0(O) is not defined, since G" L^Q; V). To complete the argument, we show that, 

for this particular kernel, the linear map A is injective. Indeed, let G A be such that 

-i-t 

■ ■ "is. 



= A0(t) = f (1-t- s)0(s)di 
Jo 



The above equality readily implies that = 0. 

Let us provide examples also for (ii) and (iii). Again, u G V is any nonzero vector. 
Example 10.6. With // as in the previous example, set 



F(t) 



1 

sin sm 1 + t cos 1 — cos 1 

t 



X[o,i](t) u. 



Then, F = A0 with 



0( S ) 



cos 



i - s) 3 i - s (i - s y 



sm 



'U. 



Note that ||F|| V G L°°(M + ) but lim t ^ F(t) does not exist in V. 
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Example 10.7. Consider the kernel 



M«) = \/— » « = (o,i), 



and set 



Then, F = A0 with 



I ft 



F (t)= I \l g (t ^ g) S ds X[o,i](*K 



0(s) = -^=U. 

V s 

It is easily verified that is summable on IR + but 

Hm||F(0||v = oo. 

Remark 10.8. A related (and very challenging) question is the following: given a function 
F G Co([t, oo), V) for every t > 0, find (easy to handle) conditions ensuring that F £ S. 
In fact, the answer seems to be strongly dependent on the particular choice of the kernel. 
For instance, with fx as in Example 10.5, F <E S if and only if F(t) = DF(t) = for t > 1, 
F G L X (]R + ; V), and DF is absolutely continuous from [t, 1] into V for all t > 0. In which 
case, F = A0 with 0(s) = _D 2 F(1 — s). On the other hand, if fx(s) = e~ s , then F G S if 
and only if F(t) = e - *« with u G V (cf. the next Example 11.7). 

11. Proper States: Recovering the Original Equation 

The purpose of this section is to establish the link between (8.1) and the original equation 
(3.1), up to now only formal. To this end, we have to recall the particular form of the 
initial datum £o> obtained in a somewhat heuristic way in Section 6. This gives a clue 
that not all the states are apt to describe the behavior of the original equation, but only 
certain particular states having a well defined structure. 

Definition 11.1. A vector £ G V is said to be a proper state if 

S(t)=DF{t\ 

for some F G S. We denote by V the normed subspace of V (with the norm inherited by 
V) of proper states. 

For any given kernel /x, an immediate example of proper state is 

f (r) = fi(r)u, ueV. 

Indeed, £ = DF with 

F(t) = -M(t)u = - fx(t + s)u ds. 
Jo 

Lemma 11.2. Let £ G V . Then, there exists a unique F G S such that £ = DF. Besides, 
F belongs to S . Moreover, for every <fi G A such that F = A(p, it follows that 



e(r)= /V(r + s)<t>(s)ds. 
Jo 
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Conversely, if £ G V has the above representation for some G A, then (G? and 

£(r) = DA<j>(r). 

Proo/. Let F E S be such that f = DF. From Lemma 7.3, L>F G L X (Q; V). Hence, the 
map 

ri 



-J DF(r)dr = F(t) 



belongs to Co([0, oo), V). Therefore, F G <So, an d it is apparent that F is uniquely 
determined by DF. The remaining assertions follow by (10.1). □ 

Here is a concrete application of the lemma. 
Example 11.3. Let ji and u as in Example 10.7, and define 



Then, F = A0 with 



<f,(8) = 



1 1- 

—7= sin - | u. 

V s s - 



Setting 



1-T 



— sm - ds 

y/s(l -T-S)(T + Sf/ 2 S . 



u, 



it is not hard to verify that £ G V. From Lemma 11.2, we conclude that £ = DF G V. 
Note that, as expected, F E So- Indeed, 



lim F(t) = hu in V, 

t->o 



with 



x = lim 



N 



1 x-l 



N^co J 1 x V £ 

In particular, Lemma 11.2 says that the map 

T : V -> 5 

defined as 



sinrrdx ~ 0.28. 



W) = - J t(r)di 



is injective. Since 



c 5b, 



and the inclusion So C S can be strict, the map V is not, in general, onto. In fact, the 
inclusion TV C <S can be strict either. 
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Example 11.4. Let 

fx(s) = l-s, ft = (0,1). 
Given any nonzero vector u G V, consider the function 

F(t) = (y/t-i) 2 Xm (t)u. 

Then, F = A0 with 



<K*) 



.2(1 -s) 3 / 2 J 
We conclude that F G S . On the other hand, 

R -1 



u. 



DF(t) = 

which does not belong to V. 

We have now all the ingredients to state the main result of the section. 
Theorem 11.5. Let (uo,Vo,F ) G V x H x S . Assume in addition that 

F G TV. 

Then, a function u is a solution to (3.1) with initial state (u ,v ,F Q ) (according to Defi- 
nition 10.1) if and only if 

(u(t),d t u(t),?) = S(t)(u ,v ,to), 

with £' as in (8.6) with 

Ut)=Kt)u + DF (t). 

Conversely, if u is a solution to (3.1) with initial state (uo,vo,F ) and F G" TV, then 
there is no corresponding solution in the extended state space. 

Proof. Since u G C([0, oo), V), arguing as in the proof of Lemma 10.3, the equality 
J J fj,'(r + s)u(t — s)ds^dr — — J fj,(s)u(t — s)ds 

holds for every t > 0. Thus, using (8.6), keeping in mind the particular form of £ an d 
the fact that F G <So, we readily get 

(11.1) / £*(r)dr = M(0)u(t) - [ fi(s)u(t - s)ds - F (t) . 

Jo Jo 

This equality, in light of (3.4), (4.1) and (8.1), proves the first statement. 

To prove the converse, assume that u(t) is at the same time a solution to (3.1) with 
initial state (u ,v , F ), and equal to the first component of S(t)(u , v , £ ), for some 
£o £ V. We reach the conclusion by showing that Fq G TV. Indeed, calling now the 
third component of S(t)(uo, vq, £o), from (3.4), (4.1) and (8.1), we obtain again (11.1). 
Since, by (8.6), 

/ £'(r)dr = / £o(r)dr + M(0)«(t) - M(t)uo - / Ai(s)u(t - s)ds, 
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we conclude that 

[fi(r)u - eo(r)]dr = M(t)u - J e (r)dr = F (f). 
Hence, 

-h(t)u + £ (t) = DF {t), 
meaning that £o _ H u o £ V and F = r(£ — f^u ). □ 

Remark 11.6. Since we have an existence and uniqueness result in the extended state 
space, Theorem 11.5 provides an existence and uniqueness result for (3.1), according to 
Definition 10.1, whenever we restrict to initial states with F E TV. 

However, there are situations where the equality S = TV holds, as in the case of the 
exponential kernel. 

Example 11.7. For a > and k > 0, consider the kernel 

^(s) = ae" KS , Q = M + . 

Since 

/j(t + s) = e~ Kt /j,(s), 

it is apparent that 

S = S = {F(t) = e~ Kt u with u e V}. 

In turn, 

V = (£(r) = e~ KT u with ueV). 

Clearly, S = TV. 

Remark 11.8. Incidentally, the above example sheds light on another important issue: 
there exist states which are not proper states; in other words, the inclusion V C V is 
strict (and not even dense). 

In summary, there might be state functions of the original approach that have no 
corresponding (proper) states. Conversely, only the proper states describe the original 
problem. In this respect, the state approach is a more general model, which is able to 
describe within the formalism of semigroups also a certain class of Volterra equations with 
nonautonomous forcing terms. 

Nonetheless, if we start from a proper state, it is reasonable to expect that the evolution 
remains confined in the space of proper states. To this end, let us define the extended 
proper state space as 

% = V x H x V, 

which is a normed subspace of 21. 

Proposition 11.9. If z E it follows that S(t)z E 2t p . 

Proof. Let z = (wo,^o,£o) £ %p- Then, £ = DF for some F E S. In turn, F = A(f> for 
some E A. Denoting as usual S(t)z = (u(t), v(t), £*), and setting 

^\s) = u(t - s)x(p,t)(s) + u X(t,e)(s) - u(t) 

and 

0*( s ) = (f)( s -t)x(t,e)(s), 
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the representation formula (8.6) can be equivalently written as 



?(t)= /V / (r + s)K(s) + 0^)}d S . 
Jo 



By Lemma 11.2, in order to prove that G V, we are left to show that ip* + 0* G A. 
Indeed, since HS^zHa < 

ri 

^(s)||^(s)|| v ds<3M(0)||z|| a . 



J" 

Jo 



Therefore, tp* G L^(Q; V) C A. Concerning 0', we have 

fi(s)(j) t (s)ds = A(f)(t) G V, 

'o 

which yields 0* G L^(f2; 7)cA □ 

In particular, from Theorem 8.1 and Theorem 9.1, we have the following corollary. 

Corollary 11.10. The restriction 

S P (t) = S(t) l%> : % -> % 

is a contraction semigroup on Assuming also condition (9.1), £/ie semigroup S p (t) is 
exponentially stable. 

One might ask whether V (and in turn 2l p ) is a Banach space. This is true, for instance, 
for the exponential kernel of Example 11.7. However, in general, the answer is negative. 

Example 11.11. Take the kernel 

H(s) = l-s, H = (0,1). 

Let £ : [0, 1] — > [0, 1] denote the famous Vitali-Cantor-Lebesgue singular function, and 
let <L n be the usual approximating sequence of absolutely continuous functions (cf. [38]). 
Consider the sequence 

F„(t)= / (t n (s)ds X[o,i](t)u, 
l Jo J 
where u G V is any nonzero vector. Then, F n = A<f) n with 

(f) n (s) = £' n (s)u. 

Setting 

f„(r) =DF n (r) = -€ n (l - r) u, 
it is readily seen that £ n G V, and, consequently, £ n G V . It also apparent that 

lim i n = i in V, 

where 

£(t) = -C(1-t) u . 
However, £ $.V. Indeed, if not so, the function 



F(t) = rf(t) = [ / €(s)ds\ X[ o,i](t)u 
l Jo J 



A NEW APPROACH TO EQUATIONS WITH MEMORY 33 

belongs to S. Hence, there is E A such that 

K<f>{t) = J (jf (f>(a)da)ds = J €(s)ds u. 
But this implies that 

<f>{s) = €{s) u = 

for almost every s E (0, 1). Thus, F = and £ = .DF = 0, leading to a contradiction. 

12. State versus History 

We finally turn to the main issue that motivated this work: the comparison between the 
past history and the state approaches. We begin to show that each element of M gives 
rise to a proper state, defining the linear map 

T1:M 



as 



rt 

Th]{r) — — /i'(r + s)r](s)ds. 
Jo 



'o 

Lemma 12.1. Let rj E M. Then, the vector 



U.rj(r) — — /i'(r + s)r)(s)ds 
Jo 

belongs to V . Moreover, 

l|n?7llv< \\v\W 

Proof. Let r] E M. Then, 

\\U V (r)\\l<{^-^(r + S )Hs)\\ v d s y 



Therefore, 



< [ e - f l , (T + S )ds I -^(T + s)Hs)\\ 2 v ds 

Jo Jo 

= Mr) [ e -^(r + s)Hs)r v ds. 
Jo 

\\m\ 2 v< fdr [ e -^(r + s)Hs)r v ds 
Jo Jo 

= [{f G -^\r + s)dr)\\r 1 {s)\\lds 



\r>\\ 2 



Thus n?7 E V, and the norm inequality stated above holds (in fact, equality for rj(s) = u, 
with u E V). Since M. C L l (Q,\ V), because of the straightforward estimate 

WvWm, 
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and LUfl; V) C A, it follows from Lemma 11.2 that Hrj is a proper state. □ 

Rephrasing the lemma, II e L(Ai, V); namely, II is a bounded linear operator from M. 
into V. Moreover 

l|n||i(^,v) = i- 

We now clarify the correspondence between rj E M and its related proper state Urj. 
Letting 

z = (u Q , v , i] ) eVJt, z = (u , v , Tlrio) G %, 

and denoting 

E(t)z = (u(t),dtu(t),ff), S p (t)z = (u(t),d t u(t),Z t ), 
we have the following result. 

Proposition 12.2. The equalities 

u{t) = u(t) and = Tiff 

hold for every t > 0. 

Proof. Introduce the function (cf. (5.4)) 

u(t) -u(t- s) < s < t, 

r)o(s — t) + u(t) — u s > t, 

which solves the Cauchy problem in M. 

(tV = T7i t + d t u(t), 

(12.1) I dt 1 1 w ' 

[?7 = i] . 

The representation formula (8.6) for furnishes 

£*(r) = ILifa^ + r) + ^(r)u(t) - fi(t + r)w + / //(r + s)u(t - s)ds = Urj^r). 

Jo 

Thus, exploiting (10.1), 

pi pi pt 

\ e(r)dr = / ny(r)dr = At/(0) = / ^(s)^( S )d S , 
Jo Jo Jo 

and, consequently, 

pi pi 

(12.2) + / pL(s)rj t (s)ds = d tt u + A u + / £'(r)dr 

Since 

■u(O) = Mo and 9 t w(0) = i>o, 
collecting (12.1)-(12.2) we conclude that 

{u{t),d t u{t),rf) = E(t)z = {u{t),d t u{t),ff)- 
This finishes the proof. □ 



0. 
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Nonetheless, in general, the map II : M. — > V is not injective. This means that two 
different initial histories may entail the same initial proper state, so leading to the same 
dynamics in the future. 

Example 12.3. Let N G N. Given a n > and k n > . . . > k\ > 0, consider the kernel 

N 



H(s) = ^a n e KnS , n 



71=1 



For x m G R to be determined later, define 



77 (s) = U, T) N (s) = J ^2 X mS m 



u, 



m=l 

where n 6 K is a fixed nonzero vector. Clearly, r)o,r)N G and 770 7^ 77at. Besides, 

n=l 



and 



having set 



with 



N 



n=l 



N „™ TV 

7. 7 , 



s m e" K "'ds = Yl b 

in 1 m=l 



mi 



The determinant of the matrix B = {b nm } is given by 



n r n (i 4)^o. 



det x 

l<n<Af l<m<n<iV 

Therefore B is nonsingular, and we can choose a; = [xx, . . . ,x N ] T to be the (unique) 
solution to the linear system 

Mx = [1,...,1] T . 

In which case, J n — 1 for all n, so that the equality II770 = II77JV holds true. 

However, for the kernel of Example 12.3, one can verify that II maps M. onto V . Thus, 
every proper state is realized by a history from M.. On the contrary, the next example 
describes a situation where the map II is injective on Ai, but ILM is strictly contained 
in V, meaning that all different histories in Ai lead to different proper states, but there 
are proper states which do not come from histories. We need first a definition and some 
preliminary results. 
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Definition 12.4. A positive sequence {«;„}, n G N, is called a Miintz sequence if K n f oo 
and 

oo 

— = 00. 

n=l Kn 

Given a function <? G Li oc ([0, 00)) such that s t— > e~ Xs g(s) G L 1 (M + ), for some A > 0, 
we denote its (real) Laplace transform by 

/•OO 

AKx) = / e- M ff (s)ds. 
Jo 

A celebrated result due to C. Miintz says that if {n n } is a Miintz sequence belonging to 
the domain of Cg and 

= 0, VnG N, 

then <? is identically zero (see [50]). 

The following lemma is standard. A three-line proof is included for the reader's conve- 
nience. 

Lemma 12.5. Let n n > be strictly increasing, and let f3 n G R be the general term of an 
absolutely convergent series. Consider the function h : [0, 00) — > R defined as 

00 

h(t) = J2(3ne- Knt . 

71=1 

Then, h is identically zero if and only if (5 n — for every n. 
Proof. One implication is trivial. If h = 0, we have the equality 

0= e KlT h(T)dT = + V ^ f 1 - e'^- K ^), Vt>0. 

The uniform boundedness of the series forces f3i = 0. Iterate the argument for all n. □ 

We are now ready to provide the aforesaid example. 
Example 12.6. Consider the kernel 

00 

fj,( S ) = J2 a nZ~ K " S , ^ = R + , 

71=1 

with K n > strictly increasing and a n > such that 

00 

n=l 

Such a pi is summable on R + . We first observe that if g G L^(R + ), then 5 G Ll oc ([0, 00)) 
and belongs to the domain of Cg. Let us extend in the obvious way the map IT to 
the domain 

M* = {v G Lj(fi; V):t~-J //(t + s)77(s)d S G v} 
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(we keep calling II such an extension). Note that M C M.+ C A, and from Lemma 11.2 
we learn that ILM* C V. Given rj E M.+ and w E V*, we consider the duality product 

g w (s) = {r)(s),w) e Ll(R + ). 

In view of (10.1), 

Ur] = Ar] = 0. 

But At] = if and only if 

oo 

J2fin{w)e- Knt = 0, Vt>0, Vio e r, 

71=1 

having set f3 n (w) = a n Cg w (n n ). Moreover, 

OO OO „oo 

^2\(3 n (w)\ <^a n I e~ KnS \g w {s)\ds = ||^ w ||i,i(r+) < oo. 

n=l n=l ^° 

Hence, from Lemma 12.5, the above equality is true if and only if 

Cg w (K n ) = 0, VnG N, \/w E V*. 

Therefore, if {/«„} is a Miintz sequence, 

Ur] = (? m = 0, Vw E V* 7i = 0. 

In which case, the map II is injective on A4±. Accordingly, to conclude that ILM is a 
proper subset of V we have to show that the inclusion M. C M± is strict. This is obtained, 
for instance, by looking at the elements 

rj( s ) = e aKlS u, 

where a E [\, 1) and u E V is any nonzero vector. The details are left to the reader. 
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